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Abstract 

Recently, we have developed a formalism to evaluate QCD loop dia¬ 
grams with a single virtual gluon using a running coupling constant at the 
vertices. This corresponds to an all-order resummation of certain terms 
(the so-called renormalon chains) in a perturbative series and provides 
a generalization of the scale-setting prescription of Brodsky, Lepage and 
Mackenzie. In its original form, the method is applicable to Green func¬ 
tions without external gluons and to euclidean correlation functions. Here 
we generalize the approach to the case of cross sections and inclusive decay 
rates, which receive both virtual and real gluon corrections. We encounter 
nonperturbative ambiguities in the resummation of the perturbative series, 
which may hinder the construction of the operator product expansion in 
the physical region. The origin of these ambiguities and their relation to 
renormalon singularities in the Borel plane is investigated by introducing 
an explicit infrared cutoff. The ratios R e + e - and R r are discussed in detail. 
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1 Introduction 


Recently, an approach has been developed to analyse the momentum flow in 
Feynman diagrams containing a single virtual gluon line [[!]]. It was initiated 
by the observation that sometimes the “typical” momenta in a loop diagram are 
different from the “natural” scale of the process. This is indicated by the fact that 
in the perturbative series for some quantities depending on a single mass scale M , 
there remain large two-loop corrections of order a 2 when one uses the “natural” 
scale M to evaluate the running coupling constant. Here ( 3 q = 11 — | n/ is the first 
coefficient of the /3-function, and rif denotes the number of light quark flavours. 
Brodsky, Lepage and Mackenzie (BLM) have argued that the appearance of such 
terms indicates an inappropriate choice of the renormalization scale, and that one 
should eliminate them by readjusting the scale in the one-loop running coupling 
constant [|2], 0- This prescription defines the BLM scale //rlm , which may be 
interpreted as the “typical” scale of virtual momenta in Feynman diagrams. One 
of the goals of Ref. |l| was to understand the relation between the BLM scale 
and the “natural” scale of a process, in particular in cases where /U B lm -C M . 

Consider the calculation of a physical (i.e. renormalization-scheme invariant 
and infrared finite) quantity S(M 2 ) at order a s in perturbation theory. The 
BLM prescription is equivalent to using the average virtuality of the gluon as 
the scale in the running coupling constant. In Ref. |IJ we have generalized this 
proposal by performing the calculation with a running coupling constant a s (—k 2 ) 
at the vertices, where k is the momentum flowing through the virtual gluon 
line. This idea is not new (see, for instance, the discussion in Ref. ||), but 
because of its computational complexity it had not been pursued previously to 
the point of practical implementation. The result of such a calculation, which 
gives the average of the running coupling constant over the virtual momenta in 
Feynman diagrams, may be written as (for simplicity we normalize the lowest- 
order contribution to unity) 

S ies (M 2 ) = 1 + jdrw(r) , (1) 

o 

where the scheme-independent function w(t) describes the distribution of virtu¬ 
alities in the loop calculation. The constant C depends on the renormalization 
scheme in such a way that the value of the coupling constant a s (re c M 2 ) is 
scheme-independent (in MS-like schemes). This implies that the product 

Ay = e ~ c A Aq CD (2) 

is scheme-independent, where Aqcd is the scale parameter in the one-loop ex¬ 
pression for the running coupling constant. We note that C = —5/3 in the MS 
scheme, C = —5/3 + 7 — In 47r in the MS scheme, and C — 0 in the so-called 
V scheme, in which the running coupling constant is defined in terms of the 
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heavy-quark potential 0. Although the V scheme is particularly convenient for 
our considerations, we shall display the dependence on C in order to keep the 
discussion general. 

Eq. (P is equivalent to an all-order resummation of certain terms in the per¬ 
turbative series for the quantity S(M 2 ). This explains the subscript “res”. Using 
the one-loop expression for the running coupling constant to relate a s (re c M 2 ) 
to a s (M 2 ), one finds that 

S ro ,(M 2 ) = 1 + £ (^P-JcnPS - 1 , (3) 

where the coefficients c n are given by the integrals 

OO 

c n = J(1 tw(t) (-C - lnr) n_1 . (4) 

o 

In other words, our approach resums all terms of order /3 q in the perturbative 
series for the quantity S(M 2 ). The resummation of such terms has also been 
considered by Beneke and Braun (see also 0), using however a different 
formalism. To understand what it corresponds to, consider a Feynman diagram 
with a single virtual gluon line. A particular set of higher-order graphs is obtained 
by inserting (n — 1) light-quark loops on the gluon propagator. The resulting 
contributions are of order In an abelian theory they are obviously related 

to the renormalization of the coupling constant. In a non-abelian theory one 
may try to incorporate the effects of gauge- and ghost-field loops by replacing 
rif by — | (3q. After this replacement, these contributions are no longer related in 
an obvious way to subclasses of diagrams; however, they are called renormalon 
chains. It is the all-order resummation of renormalon chains that is accomplished 
in (|). 

It is instructive to relate the approximation S Tes (M 2 ) to the BLM scale-setting 
prescription We find 

SUM 2 ) = 1 + N m) jl + A + | 

= S blm (M 2 ) + @1 «s(Ablm) + - - ■ . (5) 

where 


N 


I 1 BLM 



0 


exp ((In t) + C)m 2 , 


A = a 2 = (ln 2 r) — (lnr) 2 . 


( 6 ) 
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We use the symbol 


Jdr w(t) f(r) 

(nr)) = - (7) 

/ dr w(t) 
o 

for the average of a function /(r) over the distribution w(r). The one-loop 
coefficient N, the value of the coupling constant a<j(/i| LM ), and the parameter A 
are renormalization-scheme invariant. The first correction to the BLM scheme 
appears at order (3 q a^ and is related to the width (with respect to In r) of the 
distribution function. 

Let us now investigate the infrared (IR) properties of the integral representa¬ 
tion (H|). The fact that the integration extends to r = 0 indicates the appearance 
of nonperturbative effects, which arise because of the divergent behaviour of per¬ 
turbative expansions in QCD. Because of a factorial growth of the expansion 
coefficients c n , a perturbative series such as the one on the right-hand side in 
([D has a vanishing radius of convergence. It is useful to introduce a generating 
function for these coefficients, 

oo n —1 

(8) 

which can be identified with the Borel transform of the series (|3]) with respect 
to the inverse coupling constant [16!], in the limit where /3q —s► oo. The series on 
the right-hand side of (^) is usually convergent in a finite interval u E]—j,k[ 
around the origin, where j and k are called the positions of the nearest UV and 
IR renormalons (see below). Note that we have defined the function S(u) in a 
scheme-independent way. Its relation to the distribution function w(t) is [1] 

OO 

S(u) — JdTw(r)T~ u , 

0 

mo-Moo 

w(t) = - —: [ d uS(u)r u ~ 1 . (9) 

2m J 

UQ—ioO 

The choice of the parameter Uq in the inverse Mellin representation is arbitrary 
as long as it is inside the interval ] — j, k [, where the first integral is well defined. 
Outside this interval the Borel transform S(u) is defined by analytic continua¬ 
tion. The fact that the function w(t) is well-defined and has a concrete physical 
interpretation, whereas the definition of S(u) requires an analytic continuation, 
will become important later. 

The Borel transformation can be inverted to give 

- 00 , 

1 r ^ / Air'll \ 

Sbo„,(M 2 ) = 1 + Ja Jdu S(u) e- c ” exp ( - fta>(Af2) ) 
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( 10 ) 


= 1 + 


A) 


d uS(u) 


'M 


V 


\M 2 


where Ay is the scheme-independent parameter defined in If the integral 
existed, it would define the Borel sum of the partial series on the right-hand side 
of (|^). In general, however, the Borel transform S(u) contains singularities on the 
real u-axis, and the result of the integration depends on how these singularities 
are regulated. Much of the nonperturbative structure of QCD can be inferred 
from a study of the Borel transform ||-[]T7|. Its singularities on the negative 
axis arise from the large-momentum region in Feynman diagrams and are called 
ultraviolet (UV) renormalons. They are Borel summablc and pose no problem to 
performing the integral in (0)- The singularities on the positive axis arise from 
the low-momentum region in Feynman diagrams and are called IR renormalons. 
Their presence leads to an ambiguity in the evaluation of the Borel integral, 
reflecting the fact that in (|j) one is attempting to sum up a series which is not 
Borel summable. 

In our approach, IR renormalon ambiguities appear because the r-integral 
in d||) runs over the Landau pole in the running coupling constant. In order to 
regularize the integral, one may write 


/3o Qt s {re c M 2 ) 


47T 


In M 2 /Ay + In r 


= P 


(t --) +77<5(lnr-lnr L ), (11) 

V In t — m tl / 


where tl = Aiy/M 2 is the position of the Landau pole, “P” denotes the prin¬ 
ciple value, and rj is a complex parameter, which depends on the regularization 
prescription. A particularly useful regularization is to set 77 = in, corresponding 
to 

Ao a s (Te c M 2 ) _ 1 

An In r — In tl — ie 

Let us show that the resummation (|lj) together with this regularization prescrip¬ 
tion is equivalent to performing the Borel integral ([L0[) along a contour above the 
real w-axis. We note that 


Sies(M 2 ) = 1 + [ d tw(t) --- 

Po J In r — In tl — xe 

00 00 

- 1 + L /drffi(r) [da(—y"’e- 

/3„ J { Vt l J 

_ 200 . 0 
= 1 + A/ d “ s <“K#) % ""' 

0 


(13) 


where we have used the integral representation (|9|) of the Borel transform S(u) 
in terms of the distribution function w(t). Note that interchanging the order 
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of integrations when passing from the second to the third line is justified, since 
both integrals converge. Let us now assume that the Borel transform S(u) is an 
analytic function in the complex w-plane apart from singularities on the real axis. 
It is then allowed to rotate the integration contour in the last integral in (|T^) 
into a contour that runs above the real axis, and then to take the limit e —► 0. 
This proves the equivalence of our resummation with the Borel integral in (|Tbj). 
Using the same argument it follows that choosing a +ie to regulate the Landau 
pole in ( |T2| ) corresponds to performing the Borel integral below the real w-axis, 
and regulating the Landau pole with a principle value prescription is equivalent 
to taking the principle value of the Borel integral. 

The dependence on the regularization prescription leads to an intrinsic am¬ 
biguity in the definition of S res (M 2 ). Inserting ( |TT| ) into (jlj) and defining the 
renormalon ambiguity AS'ren as the coefficient of rj, we find 




tl 

A) 


w[t l ) ~ 


Wo 

A 



(14) 


where we have used the fact that Tl -C 1 to expand the distribution function: 


w{t) 


WoT k 1 + ... for t —> 0. 


(15) 


It is the asymptotic behaviour of the function w(r) for small values of r that 
determines the size of the renormalon ambiguity. The power k coincides with 
the position of the nearest IR renormalon pole in the Borel plane [IJ. Note that 
k > 0 in order for the integral in ([!]) to be IR convergent. 

In Ref. |Tj we have developed the resummation technique for QCD Green 
functions without external gluon fields and for euclidean correlation functions of 
currents. In these cases only virtual gluon corrections have to be considered. Both 
from the phenomenological and from the conceptual point of view it is interesting 
to extend the formalism to the case of cross sections and inclusive decay rates, 
which receive both virtual and real gluon corrections. This generalization is the 
subject of the present work. In particular, we shall consider the perturbative 
series for the cross section a(e + e~ —> hadrons) and for the decay rate T(r —> 
hadrons). In the MS scheme, one finds at order a 2 


v. 


m-m 


1 cr(e + e —> hadrons) 
3(E Ql) cr(e+e~ -> n + H~) 

1 T(r —> v T + hadrons) 
3 T(r —> u T e v e ) 


1 + + ( 0 . 17 ^ + 0 . 08 ) 

7T 



2 


x , a s{m 2 r ) 
71 


+ (0.57A) + 0.08) 



() 6 ) 


where Qi denote the electric charges of the quarks, s is the centre-of-mass energy, 
and quark mass effects are neglected. The terms proportional to give the main 
contribution to the two-loop coefficients. In the case of R r they are numerically 
quite important. If one uses the BLM prescription to absorb these terms into a 
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redefinition of the scale in the running coupling constant, one obtains /Ublm — 
0.71 y/s and /ig LM — 0.32 m T (in the MS scheme), respectively. A more striking 
example is provided by the parton model prediction for the semilcptonic decay 
rate T(b uev e ). The corresponding perturbative series at two-loop order is 



£(6 

r 


u e u e ) 


tree level 


1 - 2.41 


a s (ml) 


7r 


(3.22 fa + k ) 



(17) 


where the constant k is yet unknown. The corresponding BLM scale is /ijhwi — 
0.07mfe. Here and in the case of R T1 the BLM scales are so low that one may 
doubt the reliability of the perturbative expansion. It is therefore important 
to understand the origin of these low scales and to consider higher-order terms 
in the series. The fact that the contributions of order (3 0 a 2 in flTE|) give the 
dominant two-loop corrections justifies, to some extent, that we concentrate on 
the class of higher-order contributions provided by renormalon chains. Clearly, 
the resummation in ([|) does not replace an exact higher-order calculation, but 
nevertheless it provides a nontrivial all-order resummation of a gauge-invariant 
subset of corrections. As such, it seems worth while to apply this approach in the 
cases considered above. However, quantities defined in the physical region (i.e. 
the region of time-like momenta) differ in two respects from those considered in 
Ref. [|j], which were defined in the euclidean region. First, in the calculation of 
radiative corrections both virtual and real gluons contribute, and only the sum 
of their contributions is IR finite |27| . |28| | . Clearly, in such a situation one has to 
generalize the idea of performing a one-loop calculation with a running coupling 
constant a s (—k 2 ), since k 2 = 0 for real gluons. Second, the operator product 
expansion (OPE), which provides the framework for a systematic analysis of 
nonperturbative effects, can be justified in the euclidean region only, although it is 
sometimes argued that a “generalized OPE” based on quark-hadron duality holds 
in the physical region after applying a “smearing procedure” 0 • Below we shall 
discuss in detail the resummation of renormalon chain contributions for the ratios 
R e + e — and R T1 encountering some new features arising from the above-mentioned 
complications. The “linear” form of the integral representation in (|lj) is replaced 
by “non-linear” representations, in which instead of the coupling constant there 
appears a process-dependent function of the coupling constant. The choice of this 
function is not unique, leading to nonperturbative ambiguities in the definition 
of the resummed series that are not related to IR renormalons. Based on our 
results we cannot exclude that the one-to-one correspondence between the size of 
long-distance contributions and the position of renormalon singularities, which in 
the euclidean region establishes the link between the OPE and the Borel integral, 
may be lost in the physical region. 

In Sect. |]we review some results derived in Ref. [p]] for the perturbative series 
for the correlator of two vector currents in the euclidean region. We discuss the 
connection between IR renormalons and nonperturbative effects. By introducing 
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a factorization scale, which serves to separate short- and long-distance contribu¬ 
tions, we construct the OPE to order 1/(Q 2 ) 2 . In Sect. |3]we perform the analytic 
continuation to the physical region and discuss the resummation of renormalon 
chains for the cross-section ratio R e + e -. Section f| is devoted to the analogous 
discussion for the decay-rate ratio R r . In Sect. [5] we investigate the origin of the 
nonperturbative ambiguities in the definition of the resummed series. We study 
long-distance effects and their relation to the singularities in the Borel plane by 
introducing an IR cutoff to regularize the Borel integral. A numerical analysis 
of the results is presented in Sect. |6|. In Sect. |7] we summarize our results and 
point out some phenomenological implications. The relation of our approach to 
the resummation procedure proposed by Beneke and Braun [£f| is discussed in the 
Appendix. 


2 Current correlator in the euclidean region 


Consider the correlator Tl(Q 2 ) of two vector currents j M = q^q in the euclidean 
region ( Q 2 = — q 2 > 0): 

i /d 4 x e‘i- (0 | T{f(x),j-(0)} 10 ) = («V - (/V) n(Q 2 ) . (18) 


For simplicity we shall neglect the quark masses. The momentum transfer Q 2 
Ay provides the large mass scale. The derivative of Tl(Q 2 ) with respect to Q 2 is 
UV convergent. We define 

D(Q 2 ) = (19) 

In Ref. [jl| we have constructed the resummation of renormalon chains for the 
corresponding perturbative series. The result is 

D res {Q 2 ) = 1 + J d tw d (t) as(yTe ^ ^ , (20) 

o 


where the distribution function is given by 


w d (t) = 8 C F \ ( j — lnr ) r + (1 T r) L 2 (— r) + lnr ln(l + r) 


r < 1, 


w d (t) = 8C F |l + lnr + Q + ^ lnr^j ^ 

+ (1 + r) L 2 (— r' 1 ) — lnr ln(l + r^ 1 ) 


r > 1 


( 21 ) 


Here Cf = 4/3 is a colour factor, and L 2 (x) = —/ 0 x 4^1n(l — y) is the diloga¬ 
rithm function. The function Wd{t) and its first three derivatives are continuous 
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at r = 1, but higher derivatives are not. A graphical representation of the distri¬ 
bution function is shown in Fig. |l|. We find it most useful to show the product 
twd{t) as a function of lnr, since then the integrals (liTV) have a direct graph¬ 
ical interpretation. In order to associate mass scales with the r-values in the 
figure, we note that In (/i 2 /Q 2 ) = lnr + C, where /i is the scale in the running 
coupling constant. In the V scheme, the point lnr = 0 corresponds to /i 2 = Q 2 \ 
in the MS scheme, it corresponds to lnr = 5/3. We observe that the distribution 
function is rather narrow and centred around lnr ~ 1. The long arrow indicates 
the average value (lnr). In order to calculate the BLM scale and the width 
parameter A defined in (|), we compute 


(lnr) = 4C(3)-~ 0.975, 

6 

(ln 2 r) = 18 - 12C(3) ~ 3.575 , (22) 


where £(3) ~ 1.20206. This leads to 


hBLM 

VQ 7 


~ 1.628 e c/2 


A D ~ 2.625, 


MS 


0.708, 


cd — 1.620. 


(23) 


Note that both in the MS and in the V scheme, the BLM scale is close to the 
“natural” scale \fCp. 

The asymptotic behaviour of the distribution function for r —> 0 is 

wd(t) = C F |6r + (4lnr — 6) r 2 + 0(r 3 )|. (24) 


As mentioned in the introduction, the linear term signals that the nearest IR 
renormalon singularity in the Borel transform of the perturbative series is located 
at u = 2. Indeed, in the case of the function D(Q 2 ) the Borel transform is O, [14 


_ 32 C F “ (-l) fc 
= 2^ U [k 2 -(l-u) 2 f 


6 C F 

7 ) -b ■ ■ • , 

2 — u 


(25) 


where the ellipses represent terms that are regular at u — 2. The presence of IR 
renormalons leads to an ambiguity in the value of the resummed series, which 
according to (|T4|) is given by 


A D ren = 


M 


V 


6 C f 




mw ) * if i# 


A 2 v 


(26) 


The appearance of IR renormalons acts as a reminder that the result of any 
perturbative calculation in QCD is incomplete; it must be supplemented by non- 
perturbative corrections. Only the sum of all perturbative and nonperturbative 















In x 

Figure 1: The distribution function as a function of lnr. The 

long arrow indicates the average value of lnr, which determines the BLM 
scale. The short arrows show the point r = A 2 /Q 2 for A = 1 GeV and 
Q 2 = m 2 (right) and (20 GeV) 2 (left). 


contributions is unambiguous. Unlike any finite-order calculation, the repre¬ 
sentation 0 makes explicit that perturbative calculations contain long-distance 
contributions from the region of low momenta in Feynman diagrams. Moreover, 
it provides a convenient way to implement Wilson’s construction of the OPE [|18[ . 
Since the integration variable r can be interpreted as a physical scale parameter, 
one can separate the contributions from different momentum scales by splitting 
up the integral into a short- and a long-distance piece. We define 

a s {re c Q 2 ) f ^ a s (re c Q 2 ) 

— -+ j d rw D (r) 

A 2 /Q 2 0 


D V es(Q 2 ) = 1 + J d TW D (t) 


4n 


= 1 + AdW 2 ,A 2 ) + AdW 2 ,A 2 ). 


(27) 


Here A acts as a factorization (or separation) scale, which should be chosen such 
that Ay < A a JQ 1 ■ The short-distance contribution can be reliably calculated 
in perturbation theory, and it is free of renormalon ambiguities. The long-distance 
contribution must be combined with other nonperturbative corrections. Only 
the sum of all long-distance contributions is well defined. The dependence on 
the arbitrary scale A must cancel in the final result. This A-dependence can be 
controlled in perturbation theory by means of the renormalization-group equation 


A 2 Ad(Q 2 , A 2 ) = -A 2 -A- D sd (Q 2 , A 2 ) = 


a s (e c X 2 ) A 


^w d (\ 2 /Q 2 ). (28) 


dA 2 4vr 

If A 2 /Q 2 -C 1, the Q 2 -dependence of the long-distance contribution is determined 
by the asymptotic behaviour of the distribution function for small values of r and 
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coincides with the Q 2 -dependence of the renormalon ambiguity. In the present 
case, it follows that Ad(Q 2 , A 2 ) oc 1/(Q 2 ) 2 . At the same order in the OPE there 
appear nonperturbative effects parametrized by the gluon condensate 130], which 


has an UV renormalon ambiguity that compensates the ambiguity from the IR 
renormalon in the perturbative series f§|j-[ p~j| 
for the physical correlator 


Hence, to this order we may write 


D(Q ) = 1 + D sd (Q 2 , A 2 ) + Ad(<2 , A ) + 


2tt 


3 (Q 


2\2 


(a s G~) + ... 


= 1 + Ad (Q i A ) + 


2vr 


3 (Q 


2\2 


(a s G 2 ){ A) + ... , 


(29) 


where the last equation defines a scale-dependent gluon condensate, which is free 
of renormalon ambiguities. From (^8|), it follows that 

A 2 A( Q ,G 2 > ( A) = ^£ ; Q , ( e c A 2 ) A 4 . (30) 

Let us now present a numerical analysis of our results. We choose Q 2 = m 2 
as the momentum scale. In Table [l] we compare the following approximations for 
D{m 2 )\ the one- and partial^ two-loop expressions evaluated using the “natural” 
scale Q 2 in the running coupling constant, the one-loop expression evaluated using 
the BLM scale, the truncated series including the Erst correction to the BLM 
scheme (denoted by _D B lm*) given by the term proportional to A in ((^), and the 
partial resummation of the series provided by the integral over the distribution 
function regulated with the principle value prescription. We use the one-loop 
expression for the running coupling constant with rif — 3 light quark flavours. 
We take a s (m 2 ) = 0.32 in the MS scheme, corresponding to A^ 13 ~ 201 MeV. The 
corresponding value of the scale parameter in the V scheme is A y ~ 461 MeV. 
We observe that the effect of higher-order corrections is not negligible. The 
resummation of renormalon chains leads to a 3.3% increase in the value of H(m 2 ) 
with respect to the two-loop result. This effect is twice as big as the two-loop 
correction itself. Most of the higher-order corrections are taken into account if one 
includes the first correction to the BLM scheme. A nice way to compare the all¬ 
order resummation of renormalon chains with the BLM approximation is to define 
a scale /i* such that the one-loop correction evaluated at this scale reproduces the 
resummed series, i.e. D res (m 2 ) = 1 + a s (/i 2 )/7r. We End /i*/bBLM — 0.69. The 
fact that /i* < /U-blm means that higher-order corrections eflectively decrease the 
average virtuality. 

As mentioned above, the problem of IR renormalons can be avoided by in¬ 
troducing a factorization scale A to separate short- and long-distance contribu¬ 
tions. For the two values Q 2 = (20 GeV) 2 and Q 2 = m 2 , the factorization point 

lr To be consistent, we take into account only the part of the two-loop corrections proportional 
to /?o (Xg. For the function D{m the remaining two-loop correction is 0.08 (a s /-7r) 2 ~ 8 x 10 -4 . 
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Table 1: Comparison of various approximations for the euclidean correlator D(Q 2 ) 
at Q 2 = m 2 . 


D\ loop 

D'2 loop 

-Dblm 

-Dblm* 

Dyes 

AD ren 

1.102 

1.118 

1.121 

1.145 

1.151 

3.3 x 10“ 3 


t = A 2 jQ 2 for A = 1 GeV is indicated by the short arrows in Fig. |I].p| Whereas 
long-distance contributions are completely negligible for Q 2 = (20 GeV) 2 , they 
are important for Q 2 = m 2 . This is reflected in Fig. |2j, where we show the 
short-distance contribution 1 + D sd (m 2 , A 2 ), defined in O. as a function of the 
factorization scale. The solid line is used in the range Ay < A < m r , where 
D sd (m 2 ,X 2 ) is well defined. The dotted line is used for the region below the 
Landau pole. Roughly speaking, perturbation theory can be trusted down to a 
value A 0 — 0.93 GeV (corresponding to /i ~ 0.40 GeV in the MS scheme), where 
a s (e c Ag) = 1. The variation of the short-distance contribution with A provides 
an estimate of the importance of long-distance effects. In the present case, these 
effects are of order a few per cent. 



Figure 2: Short-distance contribution to D(m 2 ) as a function of the 
factorization scale. The dotted line is used in the region below the Landau 
pole, where the short-distance contribution is no longer well defined. 


2 In the MS scheme, A = 1 GeV corresponds to the rather low scale /i = e 5//6 A ~ 0.43 GeV. 

The value of the ratio A/Ay = \x/ Aqcd is scheme independent. 
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3 Analytic continuation and R e + e - 


Let us now consider the analytic continuation of the euclidean correlator n(Q 2 ) 
to the physical region, where s = —Q 2 > 0. The imaginary part of the correlator 
is related to the total cross section for the process e + e _ —> hadrons at the centre- 
of-mass energy s by 


R 


\s) = 


ale' e 


hadrons) 




a(e + e~ 

where Qi denote the electric charges of the quarks. 


12ir(E <3?) Im n <-» + «) = 3(E Qi) • 


(31) 

Our goal is to resum the 
renormalon chain contributions (i.e. terms of order /5o _1 cd?) in the perturbative 
series for the quantity R(s), and to derive a representation of the result as an 
integral over a distribution function. To this end, we integrate ([20]) with respect 
to In Q 2 and use the one-loop expression for the /5-function to obtain fT] 


4tt 2 


n (Q 2 ) - n (Q 


i 91 + 1 

Qo A) 


oo 

J d twd{t) In 
o 


a s (re c Q 2 0 ) 
a s (Te c Q 2 ) ' 


(32) 


Here Q q is an arbitrary subtraction point. Next we perform the analytic contin¬ 
uation Q 2 —> — s + ie, using that for the one-loop running coupling constant 


— Im In 

7r 


1 

a s (—re c s + ie) 


Fi(a(rs)), 


(33) 


where 


and 


F\ (a) = — arccotf—^ = — arctan(7ra) + ©(—a), 
Ti \7r aj 7r 


a(^) = 


_ Po a s (e c n 2 ) 


(34) 


(35) 


47t ln^i 2 /A y 

Combining the above results we obtain an integral representation for the re¬ 
summed perturbative series R res (s), in which there appears the same distribution 
function as in (^), but a non-linear function of the coupling constant. Using the 
notation 

W e+e -(r) = w D (r) (36) 

for the distribution function in this “non-linear” representation, we write the 
result as 

1 00 

Rrll(s) = i+fl- f dr W e+e -(r) Flairs)). (37) 

A) J 

The need for the label “1” will become clear below. A similar integral represen¬ 
tation has been suggested by Beneke and Braun [|J . The precise relation between 
their approach and ours is clarified in the Appendix. These authors suggest to 
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interpret F\ (a) as an effective coupling constant, which has the attractive fea¬ 
tures that it agrees with the usual coupling constant in the perturbative region 
0 < a <C 1, but it has a smooth behaviour in the IR region. In Fig. |3| we show the 
effective coupling constant ip (a) and the bare coupling constant a as a function 
of 1/a = ln/i 2 /Ay. 



Figure 3: Effective coupling constants Fi(a) (dash-dotted line) and 
F 2 (a) (solid line) as a function of 1/a. Both functions coincide for a > 0. 

The thin line shows the bare coupling constant a. 

Since the perturbative series for R(s ) and D(Q 2 ) agree to order a 2 , the BLM 
scales are the same in the two cases. Differences appear first at order ct^, i.e. in 
the value of the parameter A in (|5|). In the case of R(s), this parameter is no 
longer given by the square of the width ar> of the distribution function. Taking 
into account that F\(a) = a — it 2 a 3 /3 + 0(a 5 ), we find that 


Ablm 


~ 1.628 e c/2 ^ 0.708, 


A e+e - = a D~Y — -0.665. 


(38) 


Since A e + e - is negative, the leading correction to the BLM scheme effectively 
increases the value of the scale. In general, non-linear representations like fl37|) , 
for which the growth of the running coupling constant in the IR region is damped, 
lead to an improved convergence of the perturbative series. 


Let us now investigate the renormalon ambiguity of the resummed series (|37| ) . 
Expanding the function F\ (a) in powers of the coupling constant and using the 
relation 

vl n+l 7T , „ , , „ d 


7T a(rs) 


=-[ (-<5r) n a(rs); S T = tt ——, 
n\ d In r 


(39) 
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we find that 


AR ien = ^~ fdTW e + e -(T)0[-8 T \8(lnT -\ut l ) = -jrO[S T \ (rW e + e -(r)) 

PO j' Po 7 r=T L 

(40) 

where r L = Ay/s is the position of the Landau pole, and 0[8 T ] denotes the 
differential operator 


0[S T ) 


, <5r ^ 

1- - H—- =F • • • 

3! 5! + 


sin 5 t 


To proceed, we use the relations 


0[8 r ] r k 
0[5 t \ r k liTY 


sin(fcTr) , k 

kn ’ 

/ d \ " sin(kn) k 

Vdfc/ kn 


(41) 


(42) 


which imply that integer powers of r in the distribution function do not contribute 
to A R ien Given the asymptotic behaviour shown in (|24j), we see that the leading 
contribution to the renormalon ambiguity comes from the term proportional to 
r 2 lnr. The result is 

Ai?„„ = -^(^) 3 + 0(l/s 4 ), (43) 


in accordance with the fact that the nearest IR renormalon pole in the Borel 
transform corresponding to the quantity R(s) is located at u = 3 ill 


32 C F sin(7ra) (—1 ) k k 4 C F 1 

7T u(2 - u) [k 2 - (1 - u) 2 } 2 3 3 -u + '"' 


(44) 


The ellipses represent terms that are regular at u — 3. 

At this point a word of caution related to the interpretation of the integral 
in ( |37| ) is in order. In the case of the linear integral representation (jl]), one can 
interpret the integration variable r as a physical scale parameter and introduce a 
factorization point r = A 2 /M 2 to separate short- and long-distance contributions. 
If one would proceed in the same way in the present case, one would conclude that 
the long-distance contribution to R(s) is of order A 4 /s 2 . Such an interpretation 
may be misleading, however. In non-linear representations there is no physical 
significance to the integration variable; for instance, it is possible to obtain differ¬ 
ent non-linear representations by changing variables or integrating by parts. An 
example will be given in the following section. In the case of (]37|), one could try 

3 However, it is known that a renormalon singularity at u = 2 appears when one goes beyond 

the large-/?o approximation §& 
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to obtain a linear representation by repeating the steps used to derive the result 
for the renormalon ambiguity given above. If W e + e - (r) was an analytic function, 
one would conclude that (1/r) 0[5 r \ [rW e + e -{r)\ is the distribution function for 
R(s) in the linear representation. Since higher derivatives of W e + e -(r) are not 
continuous at r = 1, however, one cannot perforin these steps. The resulting 
“distribution function” would contain an infinite number of 5-function type sin¬ 
gularities at r = 1. In other words, for the physical correlation function R(s) a 
linear representation of the form (|37| ) does not exist. As we will prove in Sect. 
the same is true for all quantities defined in the physical region. 

Let us now discuss an alternative method to resum renormalon chains for 
R e + e ~. It is based on a representation of the function R(s) as an integral in the 
complex plane. Using that in the real world the discontinuities of D(—s) are 
located on the positive real s-axis, and that ImII(0) = 0, one obtains 

R{s) =47rImII(-.s + ie) = -Im / — D{-s' + ie) = — <b — D(-s') . (45) 

7r J s 2m J s 


Inserting the expression (j2(]) for the correlation function and using the fact that 

a(rs) 


a(-xrs) — — . . 1 . . 

1 + a(rs) ln(— x) 

for the one-loop running coupling constant, we find 


, , f, ^ , 1 f dx a s (—xre c s) 

.) = l + /dr» D (r)- 

0 |x|=l 


(46) 


1 

= 1 + — jdTW e + e -(T) F 2 (ci(ts)) , 


where 


J? / x _ 1 f dx a 

2 a 2m J x l + aln(— x) 
l*l=i 


= — arctan(7ra). 
7r 


(47) 


(48) 


In Fig. j| the function F 2 (a) is compared to F x (a). It is surprising that the 
resummations (|3^) and (147] ) involve the same distribution function W e + e - (r) but 
different functions of the coupling constant. Clearly, both results would have to be 
equivalent if the perturbative approximation for the function D(—s) would satisfy 
the analyticity properties of the physical correlation function. However, it does 
not. Consider, for example, the “one-loop” expression D(—s) = 1 + a s (—s)/7T. 
The cut on the positive s-axis starts at s = 0 (not s > 0), and in addition there is 
the Landau pole on the negative s-axis. The same applies for the resummed series 
in (p0|); however, in this case the Landau pole of the running coupling constant 
leads to a cut on the negative s-axis. It is not difficult to show that these two 
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defects are responsible for the difference between the two resummations. Note 
that the effective coupling constants F l (a) and F 2 (a) are identical for a > 0, 
implying that (|37|) and (|47|) correspond to the resummation of one and the same 
perturbative series. The two functions differ only for negative values of a, i.e. 
below the Landau pole. The difference between the two resummations is thus a 
nonperturbative effect. We find 

- flSW = -J 0 JdrW e+e -(r) = (/A/ + 0(1/ s 3 ), (49) 

where ry = Ay/s is the position of the Landau pole. Note that this difference is 
parametrically larger than the IR renormalon ambiguity of the resummed series 
given in (|43|) . Indeed, the effect is not related to renormalons; both resummations 
correspond to the same perturbative series and thus to the same Borel transform. 
We will come back to an analysis of this new type of ambiguity in Sect. ||. 


4 Decay-rate ratio R T 


As a second example we consider the resummation of renormalon chains for the 
perturbative series corresponding to the ratio 


R r = 


T(r —> v T + hadrons) 


T(r 




ev e ) 


= 3 T(m 2 


(50) 


For simplicity we neglect quark mass effects and small electroweak corrections, 
which are however known |3l| , Our goal is to understand in physical terms 
the low value of the BLM scale given in the introduction, ^g LM — 0.32 m T (in 
the MS scheme), and to investigate the effect of higher-order corrections by per¬ 
forming the resummation of renormalon chains. As in the previous case, R r can 
be expressed in terms of the imaginary part of the correlator Ll(—s) in the phys¬ 
ical region s > 0, or as a contour integral over the ZTfunction in the complex 
plane. We will again find that the result of the resummation depends on which 
representation one chooses. 

The first method is to relate R T to an integral over the function R(s) in- 
For fixed neutrino energy E u , the invariant mass Shad of the 
u T + hadrons is given by 


troduced in (|3T|) . 
hadronic final state in the decay r 


Shad = (Pr - Pu ) 2 = m 2 r x, 


x = i _ ^ e [o, i]. 


m T 


(51) 


In the limit where one neglects the light quark masses, and as long as one is 
interested only in terms of order in the perturbative series, the hadronic 
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dynamics is described by the function -R(shad) discussed in the previous section. 


One obtains 24, 25 


1 

T(m 2 ) =2 J dx (1 — 3x 2 + 2x 3 ) R{xm 2 ), 


(52) 


where the function 2(1 — 3x 2 + 2x 3 ) is the normalized invariant mass spectrum 
at tree level. This expression allows us to understand the low value of the BLM 
scale. We expect that 


In 


m1 


In 


(^blm ) 2 


+ (In X) , 


(53) 


where 


i 

(lnx) = 2 J dx (1 


3a; 2 + 2a; 3 ) In a; 
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(54) 


o 

is the average of hi a; over the invariant mass spectrum. Hence, the BLM scales 
corresponding to R e + e -(m 2 ) and R r should obey the relation [see (|38|)1 


6blm = e 19/24 V blm - °- 738 eC/2 m r • (55) 

We will see below that this is indeed the exact result. The fact that the BLM 
scale for R r is lower than m T simply reflects the distribution of the hadronic mass 
in the decay. 

Let us now proceed to construct the resummation of renormalon chains for R r . 
Substituting either ®) or (SIP into (P) and rescaling the integration variable 
(r —> t / x ), we immediately obtain two different representations for the resummed 
series. They are (n = 1 or 2) 


1 

= 1 + / dT F " (°( Tm h) . 


(56) 


with F n {a) as given in (|34l) and (£H), respectively. The new distribution function 
is given by 


w A t ) = 2 / T (1 “ 3x2 + ) 


It is possible to perform this integral explicitly, with the result that 

W T {r) = 166V j 4 - fj r2 - || r3 - 2L 3(-r) - 3C(3) 
^r + i + U-r)) Inr + r 2 + ^ InV 


6 


/II 0 3 2 r 3 

(- V 3r H— r 2 H- 

V 6 2 3 


6 


lnr ln(l + r) + L 2 (—t) 


(57) 


17 











for r < 1, and 


v*V(r) = i + 

n ' 1 216 48r 


17 r 2 

— t -h 2L 3 (—r x ) 

12 3 v ; 


85 

36 


1 

4r 


4 r 2 

3 r + y 


L 2 (~t x )J lnr 

In r ln(l + r~ x ) 


11 o 3 2 r 3 

-h 3r H— r 2 H- 

6 2 3 


-r 




>(58) 


for r > 1. Here L^{x) = Jq L 2 {y) denotes the trilogarithm function. In Fig. |] 
we compare the distribution functions appearing in the resummations (|37|) and 
(|56|) for R e + e - and R r . Note that the area under the two curves is the same. The 
distribution function for R r is broader than that for R e + e ~ and shifted to lower 
scales, in accordance with the physical argument presented above. 



In x 


Figure 4: The distribution functions r W e + e - (t) (solid line) and r W t (t) 
(dash-dotted line) as a function of lnr. The arrows indicate the average 
values of In r, which determine the BLM scales. 


With the function W T {r) we compute the integrals [cf. 


N = -C f = 1, 
65 


(lnr) = 4£(3) — — ~ —0.608 , 
(ln 2 r) = ^--C(3) ^4.169 


_ (59) 

72 3 V ; 

For the BLM scale and the parameter A r , which determines the leading correction 

to the BLM scheme in (|5|), we obtain 


hBLM 

m T 


~ 0.738 e c/2 ™ 0.321, 
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7T~ 


A t = (In r) — (In r)-~ 0.509. 


(60) 


The asymptotic behaviour of the distribution function for r —> 0 is 

W T (r ) = 16Cf{ 4-3C(3)- ^r+ _ ^ lnr+ ^ ln 2 r^) r 2 + 0(r 3 )|. (61) 

Using (42) we hnd that the renormalon ambiguity in the value of the resummed 
series is given by 


AT = 

L — 1 ren 


8Cp 

f3 o 


, mi 4 
Af + 3 


A ^ 3 + 0(l/m?), 


m- 


(62) 


corresponding to an IR renormalon singularity located at u — 3. This can also 
be seen from the Borel transform of the series T(m 2 ), for which we obtain 


S T (u) = 


384 C f 


sin ttu ) 


E 


(-1 ) k k 


7T u( 1 - u){2 -u )(3 - u )(4 - u) ^ [k 2 - (1 - u) 2 } 


8 C f 


(3 -uy 


+ 


+ ... 


u 


(63) 


The ellipses represent terms that are regular at u — 3. 

For later purposes it will be useful to rewrite the result (^) in another form. 
Instead of performing the ^-integration over the distribution function, one can 
perform it over the function F n (a ) and write (n — 1 or 2) 


T th m l) = 1 + -T /drSJ T (r)G„(a(rmJ)) 


(64) 


with 


and 


G n (a) = 2 Jdx (1 — 3x 2 + 2x 3 ) F n ^ 


1 + a In x 


n T (r) = w d (t ). 


(65) 


( 66 ) 


Once again we have introduced a new notation for the distribution function, since 
the form of the integral representation has changed. Eqs. (|56|) and (^) provide an 
example of the fact that the integration variable r has no significance as a physical 
scale if the integral representation is non-linear. Whereas the contribution from 
the region r < A 2 /m 2 scales like A 2 /m 2 in the case of (|56|), it scales like A 4 /m^ 


in the case of (34). Nevertheless, both representations are equivalent and lead to 
the same numerical results. 

We will now discuss a third resummation method, which uses the representa¬ 
tion of T(m 2 ) as a contour integral in the complex plane. Since in the real world 
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Figure 5: Effective coupling constants G\(a) (dash-dotted line), G 2 (a) 
(dotted line) and G 3 (a) (solid line) as a function of (a) the inverse cou¬ 
pling constant 1/a, and (b) the coupling constant a. For small positive 
values of a, the three functions have identical Taylor expansions. 


the discontinuities of the function D(—s ) lie on the positive real s-axis and start 


at a value so > 0, one can integrate by parts in fl52|) to obtain 

1 [' d X 

Tim 2 ) = — Im / — (1 — 2x + 2x 3 — x 4 ) D(—xm 2 T + ie) 
7T J X 

0 

- (l — (1 — 2x + 2x 3 — x 4 ) D(—xm 2 ). 

'I J x 


1 

2iri 


1 * 1=1 


Inserting the expression (0i for the .D-function we find 


(67) 


^circle / 


1 

(m 2 r ) = 1 + — J drO T (r) G 3 (a(rm 2 T )), 


( 68 ) 
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where 


(73(a) = —% i — (1 — 2x + 2x 3 — x 4 ) 


2ni 
8 a 

7T 


1 * 1=1 


X 


1 + a ln(— x) 


sin 2p> — cup cos 2p> 


1 + a 2 ip 2 


2 V 3 
Sill (f cos — 


( 69 ) 


is one of the functions encountered in Ref. |53]. Eqs. (|64|) and (|68p provide three 
resummations of the perturbative series for T(m^) that have the same distribu¬ 
tion function O t (t), but different effective coupling constants G n (a). These three 
functions are shown in Fig. For small positive values of a they have identical 
Taylor expansions, meaning that all three resummations correspond to one and 
the same perturbative series. Outside the radius of convergence the functions 
G n (a ) are quite different, however, leading to differences in the value of the re¬ 
summed series which have the form of nonperturbative power corrections. Again 
these nonperturbative effects have nothing to do with 1R renormalons. It is not 
difficult to trace how the differences arise from the wrong analytic properties of 
the perturbative expression for the function D(—s), which has Landau pole sin¬ 
gularities on the negative s-axis and a cut on the positive s-axis starting at s = 0 
(not s > 0). We find that 


(72(0) — Gi(a) 
Gs(a) — (*2(0) 


I —2e -1 /“ + 2e~^ a - e~^ a ; a > 0 , 
\ — 1 ; a < 0 , 

4(e" 1/ “ - e~ 3/ “) 0(o). 


(70) 


For the differences between the three resummations we obtain, using (^U|) and 

m, 


T ves( m2 r) 




T res( m l) = 


TL 


-y dr W T (r) = 






jJdrV. T (r) 


TL 

32 C F 

A 


1 A?. 


4-3C(3) —y + 0{l/m^) , 


J mi 


(71) 


where tl = A‘y/m 2 is the position of the Landau pole. As in the case of R e + e -, 
the differences are parametrically larger than the renormalon ambiguity given 
in (|52|). In the present case they are of order 1/m^ and thus numerically quite 
significant; we note that the terms shown on the right-hand side in ([71]) equal 
—6.3% and 12.6%. A new feature of the present case is that the effective coupling 
constants G n (a) differ not only for a < 0 , but also in the region above the Landau 
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pole (although by exponentially small terms). This is what makes the various 
resummations differ by terms of order \jvr? T . We will argue in the next section 
that such differences can be avoided. There is a physical argument why one 
should prefer one form of the resummation, T™ c]e , over the others. However, as 
in the case of R e + e - there will remain ambiguities related to the treatment of the 
region below the Landau pole (for a < 0), which are of order 1 jm\. These effects 
are still larger than the IR renormalon ambiguity. 


5 Cutoff regularization of the Borel integral 


The goal of this section is to understand better the origin of the ambiguities 
encountered in the resummation of renormalon chains for R e + e - and R r . We 
recall that a surprising feature of these nonperturbative effects is that they are 
not related to IR renormalon singularities in the Borel plane; the differences 
between the various resummations are parametrically larger than the ambiguity 
resulting from the nearest IR renormalon. Therefore, it is tempting to conjecture 
that these differences reflect a failure of our resummation procedure, which could 
be avoided by defining the resummation in terms of the Borel integral (Jb]), in 
which case only the usual IR renormalon ambiguities appear. This conjecture may 
be wrong, however. It has been noted already by ’t Hooft that most likely the 
Borel integral is ill-defined because of strong singularities of the Borel transform 
at infinity JSj]. The presence of resonances and multi-particle thresholds in the 
physical region is in conflict with Borel summation. It is argued in Ref. || that, 
as a consequence, the Borel transform S(u ) must diverge at u —> oo stronger 
than any exponential of u. Discussions of the implications of this observation 
can be found in Refs. |34|, [35|]. We note that the strong singular behaviour for 
large values of u is not seen if the Borel transform is calculated using dimensional 
regularization. The situation is similar to power-divergent Feynman integrals of 
the type / d d k ( k 2 ) n , which vanish by definition in dimensional regularization. 
Let us therefore study the singularities of the Borel transform in a more careful 
way. 


5.1 Cutoff regularization in the euclidean region 

The proper way to analyze long-distance effects is by introducing an IR cutoff 
A, which eliminates the low-momentum contributions from perturbative calcula¬ 
tions. Let us investigate the connection between this procedure and the Borel 
integral in the euclidean region, where the OPE provides a consistent framework 
to separate short- and long-distance contributions [|l8j. As shown in (|9]). the 
Borel transform S(u ) can be defined in terms of an integral over a distribution 
function w(r), which has a concrete physical interpretation. This integral has a 
finite radius of convergence; it is well defined if u G ]— j, k [, where j and k denote 
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the positions of the nearest UV and IR renormalon singularities. The standard 
way to proceed is to define S(u) outside this interval by analytic continuation, 
in which case the divergent behaviour of the integral for large or small values of 
r is reflected in the appearance of renormalon singularities. The Borel integral 
([10]) has an ambiguity of order (A y/M 2 ) k due to the nearest IR renormalon. It is 
usually argued that the appearance of renormalon ambiguities signals the size of 
nonperturbative corrections, which have to be added to perturbative calculations 
in QCD [RJ]— []T2j] . In order to understand the link between the Borel integral and 


the OPE, let us define a regularized Borel transform as 


S Teg (u, To) = J d T w{t) T~ 

TO 


To = 


M 2 


(72) 


The dependence of To on M is dictated by the fact that, according to ([[]), the 
product y/rM has the interpretation of a physical scale. The introduction of this 
cutoff eliminates all IR renormalon singularities .0 In particular, if the asymptotic 
behaviour of the distribution function for r —> 0 is as shown in the behaviour 
of the regularized Borel transform in the vicinity of u = k is 


&*(«, To) = (To 1 "” - l) 


u — k 


M 2 

+ ... = w 0 In —— + 


(73) 


instead of S(u) = —w 0 /(u — k) + ..., which is obtained without using a cutoff. 
The regularization eliminates all singularities for u > 0, but at the same time it 
leads to a strong exponential growth of the Borel transform for u —> oo: 


S reg (u, t 0 ) ~ exp (u In — ) = exp ( u In 


To 


M 2 
A* - 


(74) 


This growth becomes arbitrarily strong as one tries to eliminate the IR cutoff, in 
accordance with’t Hooft’s conjecture. 

The regularized form of the Borel integral reads 


OO OO 


oo >2 oo oo 2 _ 

SS.,(M 2 ,r 0 ) = 1+j /d«S I8t («,r„) (j0 = 1+- J du J dr w(t) (^-) " 

0 0 T 0 17 

(75) 

The integral over u is convergent if To > Ay/M 2 , or equivalently if A > Ay, which 
we have to assume for any sensible IR regulator. It is then allowed to interchange 
the order of integration and perform the integral over u first, which gives 


d u 


tM 2 ' 


( M 


V 


A) a s (re c M 2 ) 


1ii(tM 2 /A 2 


v) 


4(/r 


( 76 ) 


4 A similar regularization can be applied to eliminate the UV renormalon singularities. 
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This leads to 


SSUM\t 0 ) = 1+ J&T<Hr) a ‘ (Te ^ m = l + S^(M\x\ i _ r ^, (77) 


7"0 


which is nothing but the short-distance contribution to resummed series defined 
in (^7]). The standard form of the Borel integral is recovered if one sets To = 0 and 
defines the right-hand side of ( |76|) to be the analytic continuation of the u- integral 
in the region below the Landau pole, i.e. for r < Ay/M 2 . Then the integration 
over the Landau pole leads to an ambiguity of order (A y/M 2 ) k in (|77|) , which is 
the usual IR renormalon ambiguity. However, the more physical way to proceed 
is to keep the IR cutoff. In the regularized version of the Borel integral it is the 
way in which the factorization scale appears that determines the size of long¬ 
distance contributions, which have to be added to the above result in order to 
cancel the dependence on A. From the renormalization-group analysis in Sect, 
it follows that this long-distance contribution has the same M-dependence as 
the renormalon ambiguity [compare ( fZE| ) with (^3])]. Hence, the IR behaviour in 
the cutoff version of the Borel integral is in one-to-one correspondence with the 
appearance of renormalon singularities in the Borel plane. It is this fact which, a 
posteriori, gives a physical significance to these ambiguities, and which establishes 
the connection between the Borel integral and the OPE in the euclidean region. 

Let us now investigate what changes if the quantity of interest is defined 
in the physical region. The crucial difference is that in this case the integral 
representation of the Borel transform shown in (j0|) does not exist. We will now 
prove this statement, which was already mentioned in the discussion of R e + e - in 
Sect. |j. Consider the following alternative integral representation for the Borel 
transform S(u ) of a quantity S(M 2 ): 


d sm itu r _ 

S(u ) =- / d tW (t) r 

7 tu J 


(78) 


Beneke and Braun have shown that this representation always exists [j|]. Let us 
assume that also the representation (|9j) exists. Then there is a simple relation 
between the functions w{r) and W{t). Rewriting the sin-function in terms of 
T-functions, we find that 

OO OO 

J d tW{t) t~ u = T(1 + u ) T(1 — u) J d tw(t) t~ u 


cxj i 

= J d tw(t) t~ u J dxx u {l — x)~ u , 


(79) 
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and introducing new variables z = x/( 1 — x) and t — r/z we obtain 


OO OO OO 

J drW (r) t~ u — j dt J dz “ ui(zt) t 


(i + z y 

0 o o v ' 

From this equation we read off the relation (a = zt ) 


( 80 ) 


W(t) — / dcr 


aw(cr) 


cr + r 


= / dcr-— 

J a + t da L 


a w(a ) 


(81) 


which establishes the connection between the two functions.^ Provided a function 
w(t) satisfying the integral representation ([|) exists, it follows that the function 
W (r) is analytic in the complex r-plane with a branch cut along the negative axis. 
We can now use the observation of Beneke and Braun that the function W (r) can 
be calculated by performing loop integrals with a finite gluon mass m 2 = tM 2 H. 
In the euclidean region these integrals have discontinuities only for unphysical 
values m 2 < 0, and the function W{r ) has indeed the analyticity properties 


following from (|3T1). However, in the physical region there is the possibility of 


real gluon emission if m 2 < M 2 , and hence the function W(t) must contain a 
contribution proportional to 0(1 — r), which is in conflict with (pT|) . Therefore, a 
function w{t ) satisfying the above integral relation does not exist in the physical 
region. 

It follows then that a “linear” integral representation of the type ([I]) does not 
exist, too. One cannot avoid to use “non-linear” representations, such as the ones 
encountered in Sects. |3] and [|. Because the form of the non-linear representation 
is not unique, however, it is no longer obvious how to introduce a factorization 
scale in order to separate short- and long-distance contributions. A possible way 
to proceed is to define the resummation of renormalon chains in terms of the 
Borcl integral. However, as in the euclidean region one should check that the 
resulting ambiguities from IR renormalons are in one-to-one correspondence with 
the size of long-distance contributions. To this end one should introduce an IR 
cutoff. A possibility which seems well motivated to us is first to regulate the Borel 
transform in the euclidean region, and then to perform the analytic continuation 
to the physical region. We shall now discuss this proposal for the quantities R(s ) 
and T(m 2 ). 


5 In the last step we have assumed that the quantity S(M 2 ) is at most logarithmically UV 
divergent, so that the distribution function w(t) vanishes for r = oo. 
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5.2 Cutoff regularization for R(s ) 

We start with the regularized form of the Borel integral for the correlator D(Q 2 ) 
in the euclidean region, which reads 

oo °o 2 _ 

D Bovel {Q 2 : T o) = 1 + 7^ / d « / &TW D (t) ■ (82) 

The integral over u is convergent as long as r 0 > A y/Q 2 . Let us now consider 
the analytic continuation of the above result in the complex plane, which is well 
defined if r 0 > Ay/\Q 2 \. Outside a circle of radius A y/r 0 around the origin, 
the regularized function D^ rel (Q 2 , To) is analytic in the complex Q 2 -plane with 
a branch cut along the negative axis. Using the above expression, we find that 
both methods of calculating the function R(s) - taking the imaginary part of 
fl(—s + fe), or integrating D(Q 2 ) around a circle of radius s in the complex plane 
- lead to the same result fl 

OO , x OO 

„ x If, sin 7 tu) _ . x frs\~ u .. 

R,, g ( S ,T 0 )=l + -Jdu-^jdTW D (T)^) . ( 83 ) 

0 TO v 

As long as tq > Ay/s, both integrals are convergent and we can perform the 
integral over u first. This gives 



smf7ru 
7 TU 




sin(7ru) 

ITU 


e - U /a 


F(a), 


(84) 


where a 
and 


l/ln(rs/Ay) agrees with the coupling constant a(rs) defined in fl35P , 


F(a) = — arccot ( — ) = — arctan(7ra); for a > 0 . (85) 

7 T \ttclJ tt 

Note that a > 0 if r > A y/ s. Setting now r 0 = A 2 /s with A > Ay, and writing 
the final result as a function of s and A 2 (instead of To), we obtain 


R reg {s, A 2 ) 


1 



f dr W e + e -(r) F(a(rs)). 

A 2/s 


( 86 ) 


For a > 0, the effective coupling constant F(a) coincides with the functions F\ (a) 
and F 2 (a ) in (^) and (|4^). However, if one wants to remove the cutoff one has 
to define F(a) in the region a < 0, i.e. below the Landau pole, where the integral 

6 Note that this is equivalent to using a representation of the form (j^) for the Borel transform 
of R(s), provided we identify W e + e -{r) = Wd(F) [see (^)[. This relation is in accordance with 
the fact that the Borel transform of R{s) in ( pd| ) differs from the Borel transform of D{Q 2 ) in 
(|25|) by a factor sin(7ru)/(7ru). 
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in ( |34D is ill defined. Both F\(a ) and F 2 (a) are “natural” continuations of F(a); 
one is continuous as a function of 1/a, the other is continuous as a function of 
a. The fact that the choice of the continuation below the Landau pole is not 
unique leads to an ambiguity in the definition of the resummed series. In the 
case of the physical correlator R(s), we have seen in Sect. [| that this ambiguity 
is parametrically larger than the IR renormalon ambiguity. 

In this context it is interesting to note that neither R)H{s) nor R)H(s) coincide 
with the principle value of the Borel integral. This is different from the situation 
encountered in the euclidean region, where taking the limit A —> 0 in ([77]) one 
reproduces the value of the Borel integral. In the present case, however, Beneke 
and Braun have shown that the principle value of the Borel integral is given by 
(see the Appendix) [ff[] 


1 1 

RBovei(s) = 1 + -p- j dr W e + e -(r) Flairs)) + — Re J dr W e + e -(r — ie) , (87) 


A) 


—TL 


where ry = A y/s is the position of the Landau pole. The imaginary part of the 
second integral determines the renormalon ambiguity: 


A R rm = 


7 T /? 0 


u 

Im J dr W e + e - (r — ie) 


-tl 


Comparing ([87]) with (^7|) and (|47|), we find that 


( 88 ) 


1 r 8CC /A 2 \ 2 

RBorei(s) - Rill(s) = -Re j &TW e+e -{T-ie) = —f{-f) +°( 1 /« 3 ), 

-TL 


1 

Po 


tl 


RBorei(s) - R[H(s) = A- Re / dr W B + e -(r - ie) 


-tl 


8 C f 




3/do 


Moreover, we obtain 


A R rm = 


4 Cf 

(3 o 


{tl - ln(l - Tl )j , 


(89) 


(90) 


which, when expanded in powers of Tl = A y/s, reproduces (|43|) . The ambigu¬ 
ities arising from the choice of the continuation below the Landau pole lead to 
differences of order 1/s 2 , which are parametrically larger than the IR renormalon 
ambiguity. 

It is not difficult to understand the size of these effects by considering the 
regularized form of the Borel integral given in (P5|). In fact, instead of taking the 
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limit A —> 0, the more physical way to proceed is to keep the IR cutoff in the 
perturbative calculation, and to add to the regularized result a scale-dependent 
nonperturbative contribution, so that the dependence on the cutoff cancels in the 
sum. From the renormalization-group equation 

A 2 A fi [8g ( s , A 2 ) = -i-F( tt (A 2 ))/V e+e -(A 2 / S ) 

= -^F(a(A 2 ))(L/ + o(l/ S 3 ), (91) 

it follows that this nonperturbative contribution must scale like 1/s 2 . At this 
point we have to add a word of caution, however. Contrary to the euclidean re¬ 
gion, one has to be careful when interpreting A as a physical scale parameter. Al¬ 
though the introduction of the IR cutoff certainly removes all soft contributions, 
we cannot exclude that it actually removes too much, i.e. contributions which 
may come from short distances as well. Unless one is able to construct a cutoff 
regularization for which the nonperturbative contribution scales like 1/s 3 , how¬ 
ever, one should be conservative and assume that there are indeed long-distance 
contributions of order 1/s 2 . 


5.3 Cutoff regularization for T{m 2 ) 

Let us now turn to the more complicated case of the function T(m 2 ). If we work 
with the representation of this quantity as an integral over R(s) as shown in (j5^) 
and insert the regularized expression given in ([861), we find 


± . , v CAJ O 

m , 9 , 9 . 9 f, sm(7ra) f , . . . fTxmz \~ l 

T reg (m;, A 2 ) = 1+ — J dx(l-3x 2 +2x 3 ) Jdu —^— j d rw D (r) y ^2 ) 


\ 2 /xm 2 


(92) 

Note that in order for the ^-integral to be convergent we have to introduce an 
^-dependent cutoff r 0 = A 2 /(xml) = A 2 /sh a d with A > Ay, where Shad is the 
invariant mass of the hadronic final state. It is thus convenient to introduce a 
new variable t' = tx and to perform the integral over x at fixed t' to obtain 


Tre g (m r ,A ) 


. 1 f , sin(7 tu) r . . TTr . .. 

l+ jJ Au ^r J At WAt \k 

0 A 2 /m? 


r'm 2 N ~ l 


v 


UXJ 

= 1 + ^ J d TW r (T)F(a(Tm 2 T )), 


(93) 


A 2 /m. 2 


where W T (r) is the distribution function defined in ( |57| ), and F(a) has been 
given in (|85|) . The freedom in the choice of the continuation to the region below 








the Landau pole (where a < 0) introduces an ambiguity of order 1 /m 2 , since 
W t (t ) —> const, for r —> 0. Two examples of such continuations are provided 
by the quantities T^J (rn 2 ) and T/A (rn 2 ) in (|56|), which indeed differ by terms of 
order 1 /rn 2 r [see ([711) 1 . If the IR cutoff is kept, the ambiguity of order 1/m 2 is 
reflected in the fact that 


1 A 2 „,, l2 


A dA^ T - g (m 2 , A 2 ) = -- F(a(X )) — W T (X 2 /m 2 ) 


A) 
16CV 
A) 


mi 


A 


4 —3C(3) F(a(X))—- + 0(l/m*), (94) 


mi 


which means that one has to add a nonperturbative contribution of order 1/m 2 
to cancel the dependence on the cutoff. 

Let us now consider the representation of T(m 2 ) as a contour integral given 
(0). Inserting there the regularized expression for the function D(Q 2 ) from 


m 


fl82[), we find 


OO OO 

T“ cle (m 2 ,A 2 ) = l+j^~ / ^ (1-2x+2x 3 -x 4 ) J du Jdrw D (r) 

° |ar|=l 0 to 

(95) 

where To = A 2 /|icm 2 | = A 2 /m 2 . This means that the integral over x can be 
performed without affecting the value of the IR cutoff, which is determined by 
the radius of the contour in the complex plane. The result isQ 


rxm T ^ u 


pcircle 

reg 


1 

K ' A2) = 1+ a/ 


du 


sin [iru 


12 


7 TU (1 — u) (3 — u) (4 — u) 


dr wd(t) 


A 2 /m 2 


rm; 

aT 


2 \ —t 


(96) 

As long as A > A y, the integral over u is convergent and can be performed first. 
Expressing the result in terms of an effective coupling constant G(a) as in (|64|), 
we find that 


T 


circle 

reg 


(m 2 , A 2 ) = 1 + 


A) 


dr fl T (r) G(a(rm 2 )), 


A 2 /m 2 


where f) r (r) = wjj(t), and 


(97) 


00 / \ 

G( a) = f du - 12 -- e-«/“ = G 3 (o) ; for a > 0, (98) 

J nu (1 — U){6 — u)(4 — U) 

7 The same result is obtained if one performs first the r-integral in (p2|), which is convergent 
for zi < 1, and then analytically continues the result to arbitrary values of u. However, this 
procedure involves an integration over regions where the Borel integral is ill defined. 
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with 63 (a) as defined in ([L]). Now the freedom in the choice of the continuation 
to the region below the Landau pole introduces an ambiguity of order 1 jm\, 
since fl T ( r ) oc r for r<l. An example of such a continuation is provided by the 
expressions for T^ cle (m 2 ) given in (| 68 |). If the IR cutoff is kept, the ambiguity is 
reflected in the fact that 


A' 


A T ™l. (m 2 jA 2 ) = _1_ G(a(X 2 )) A & T (\ 2 /ml) 


Vo 
66V 
Vo 


rnt 


G(a(A 2 ))(^)Vo(i/mJ), (99) 

\ 11 It / 


which means that one has to add a nonperturbative contribution of order lfm\ 
to cancel the dependence on the cutoff. 

It is again instructive to compare our result for T™ cle (m 2 ) with the principle 
value of the Borel integral, which according to Refs. 0,0 can be written in the 
form 


1 ^ 

T Boiel (ml) = 1 + — J dr W t (t) Fi(a(rm 2 )) + — Re jdrW T (T-ie ), ( 100 ) 


with t l = hy/m 2 r . Using (fTI|), we find that 


~ t L 


1 T L . 00 3 

T B „ r =,K) - ^ Re j AtW t (t - ie) - - j dr Si T (r) (^ - I|) 


-t~l 


48 C F /A? x 2 


Vo vm? 




Vo 

' 1/ 

+ 0(l/m®) ~ 0.03 . 


( 101 ) 


Although of order 1/m/., this difference is sizeable because of the large numerical 
coefficient. For completeness, we also quote the exact result for the renormalon 
ambiguity. It is 


4 C F 


Vo 


47 


AT ren = —— \ t l + -t l - — - AL 2 (t l ) - 2 r L - 


6 3 

13 5 


,3 T L 


In t l 


+ (i + y^n -3 ^l + y)( 1 -^) ln(l — r n)|, ( 102 ) 


which for r^Cl reproduces fl62|). 

The discussion of this section illustrates that the choice of implementing the 
IR regulator requires some amount of ingenuity. In the first case, by introducing 
the regulator we have excluded a contribution of order 1/m 2 from the pertur¬ 
bative calculation of Tim 2 /), whereas in the second case we have found a better 
regularization scheme, in which this contributions is of order 1/m/. It is impor¬ 
tant that there is a physical argument which favours the second method over the 
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first one. The integration along the circle in the complex plane involves larger 
momenta than the integration along the cut, and thus less is eliminated by in¬ 
troducing the IR cutoff. Clearly, the question arises whether there exists an even 
better representation, in which the excluded contribution is of order 1/m® and 
thus of the same order as the renormalon ambiguity. Although we cannot give 
a definite answer to this question at present, we doubt that such a representa¬ 
tion exists; we find it suggestive that a linear integral representation of the form 
(®, which would have this property, does not exist in the physical region. As 
long as one does not succeed in constructing such a regularized representation 
of the Borel integral, one has to take seriously the possibility of nonperturbative 
corrections to R T of order 1/m*. 

To summarize this section, we repeat that the new kind of nonperturbative 
ambiguities encountered in our analysis is related to the singular behaviour of 
the Borel integral for u —>■ oo. We have argued that this behaviour should be 
regulated by introducing an IR cutoff. To the perturbative calculation one then 
has to add a nonperturbative contribution in order to cancel the dependence on 
the cutoff in the final result. For quantities defined in the euclidean region, the 
power behaviour of the nonperturbative contribution is in one-to-one correspon¬ 
dence with the size of the ambiguities due to IR renormalons, which appear if 
the Borel integral is evaluated in dimensional regularization without using an 
IR cutoff. This implies that in these cases the singularities at infinity do not 
spoil the applicability of the OPE. For the cases of correlation functions in the 
physical region, however, we find that the singular behaviour at infinity leads to 
ambiguities that may be larger than those indicated by the positions of IR renor¬ 
malon singularities. At least, we could not find a cutoff regularization for which 
these ambiguities are of the same order as those arising from renormalons. If this 
conjecture is correct, it would imply that the OPE does not provide a consistent 
framework for the analysis of nonperturbative corrections in the physical region. 


6 Numerical analysis 

We now turn to the numerical analysis of our results. In Tables |2] and |3] we show 
various approximations for the quantities R{m^) and T’(m^). Compare first the 
results for R(m j:) with those for D(ntf.) given in Table [l]. The corresponding 
perturbative series are the same up to order a^. In the case of the euclidean 
correlation function, the resummation of renormalon chains leads to an increase 
of 3% with respect to the two-loop result, whereas the value of the physical 
correlation function is decreased by 1%. This reflects the smooth behaviour of 
the effective coupling constants F n (a) in the IR region. As a consequence, after 
the resummation the result lies significantly below the BLM prediction. Note 
that the difference between the two resummations R^H and R$ s is numerically 
small, although larger than the renormalon ambiguity A R ren . We also quote the 
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Table 2: Comparison of various approximations for the quantity R(s) at s = rn 2 . 


Rl loop 

R2 loop 

-Rblm 

-Rblm* 

rA) 

°res 

i? (2) 

^res 

7?Borel 

A-Rren 

1.102 

1.118 

1.145 

1.115 

1.107 

1.105 

1.105 

-6.2 x 10" 5 


Table 3: Comparison of various approximations for the quantity T(m 2 ). 


T\ loop 

T r 2 loop 

Tblm 

Lblm* 

7"(i) 

res 

jA 2 ) 

res 

o^circle 

res 

-^Borel 

AT 

— i L ren 

1.102 

1.155 

1.213 

1.238 

1.151 

1.097 

1.199 

1.228 

1.4 X 10“ 3 


principle value of the Borel integral, i?Borei, which according to (|89f ) is very close 


to RW. 

Consider now the results for T(m 2 ), which has a larger two-loop coefficient 
than D/m 2 ) and R(m 2 ). The resummation of renormalon chains again leads 
to values that are lower than the BLM approximation; however, as discussed 
in the preceding sections there are significant differences between the various 
resummations. Some of these differences are of order 1/m 2 and thus much larger 
than the renormalon ambiguity. We have also seen that, for a physical reason, 
the resummation T™ cle is to be preferred over the resummations Tj/J and Tj 2 /. 
On the other hand, the difference between T™ cle and the principle value of the 
Borel integral constitutes a real ambiguity, which cannot be avoided on physical 
grounds; at least at present we see no argument to prefer one result over the 
other. This difference is an effect of order l/m£, which is still larger than the 
renormalon ambiguity. 

Let us finally study the regularized integrals R ieg (m 2 , A 2 ) and T reg (m 2 , A 2 ) 
in (|$6D and (|97j) as a function of the cutoff scale. The corresponding curves are 
shown in Figs. |S] and jTj. The regularized integrals are well defined down to a value 
A = Ay ~ 0.46 GeV, where one reaches the Landau pole in the running coupling 
constant. Below this value we define the integrals by using one of the effective 
coupling constants encountered in the preceding sections (Fi(a) and F 2 (a) in the 
case of R(m 2 ), and G$(a) in the case of T(m 2 )), so that for A = 0 we recover the 
results given in Tables |2] and |3j. In each figure the arrow indicates the principle 
value of the Borel integral. The variation of the results in the low-momentum 
region, as well as the difference between the resummations and the principle value 
of the Borel integral, provide an estimate of the importance of nonperturbative 
contributions. 
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Figure 6: The regularized integral R, reg (m 2 , A 2 ) as a function of the cutoff 
scale. For values A < Ay — 0.46 GeV we define the integral by using the 
effective coupling constants F { (a) (upper curve) and F 2 (a) (lower curve). 

The arrow indicates the principle value of the Borel integral. 

7 Conclusions 

Generalizing the approach proposed in Ref. |T|], we have investigated the resum¬ 
mation of renormalon chain contributions in the perturbative calculation of cross 
sections and inclusive decay rates, and in particular for the ratios R e + e - and R r . 
Our approach is equivalent to an all-order resummation of the terms proportional 
to /3o _1 a;” in the perturbative series for these quantities. It provides a general¬ 
ization of the BLM scale-setting prescription [2| |3|] . in which terms of order /3 0 ct 2 
are absorbed by a redefinition of the scale in the running coupling constant. 

The discussion in Ref. [I] was devoted to QCD Green functions without exter¬ 
nal gluon fields and to euclidean correlation functions of currents, which receive 
only virtual gluon corrections. In the calculation of cross sections and inclusive 
decay rates, on the other hand, both virtual and real gluons have to be consid¬ 
ered. As a consequence, in the expression for the resummed series, which has 
the form of a one-dimensional integral over the running coupling constant with 
some distribution function, there appears an effective coupling constant, which 
is screened in the low-momentum region. Such non-linear representations have 
also been derived, using a different formalism, by Beneke and Braun Q. 

By summing an infinite set of diagrams our scheme reaches beyond perturba¬ 
tion theory. In any finite-order perturbative calculation nonperturbative effects 
are implicitly present due to low-momentum contributions in Feynman diagrams, 
but are not visible as they are exponentially small in the coupling constant. Yet 
perturbation theory “knows” about these contributions in the form of IR renor- 
malons, which make a perturbative series non Borel summable. This means that 
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Figure 7: The regularized integral T r c e g Cle (m 2 , A 2 ) as a function of the 
cutoff scale. For values A < Ay ~ 0.46 GeV we define the integral by 
using the effective coupling constant G$(a). The arrow indicates the 
principle value of the Borel integral. 


attempts to resum the series will lead to unavoidable ambiguities. In our scheme 
these ambiguities arise from the integration over the Landau pole in the running 
coupling constant. They are a reminder of the fact that perturbation theory is 
incomplete; any perturbative calculation must be supplemented by nonperturba- 
tive contributions. For the quantities considered in Ref. [jlj, the OPE provides 
the framework for a consistent incorporation of nonperturbative contributions. 
In this case the IR renormalon ambiguities can be absorbed into the definition 
of other nonperturbative parameters, such as the vacuum condensates. A new 
feature of the present analysis is that it refers to quantities defined in the physical 
region (i.e. for time-like momenta), for which the applicability of the OPE is at 
least questionable. We have found a new source of nonperturbative ambiguities 
in the definition of the resummed series, which are not related to IR renormalons. 
These ambiguities are parametrically larger than the IR renormalon ambiguities 
and thus cannot be dealt with in the standard way. 

The appearance of this new source of ambiguity is surprising and is the most 
striking observation of our analysis. In order to understand the origin and size 
of this effect, we have regularized the singular behaviour of the Borel integral by 
introducing an IR cutoff. We find that in the euclidean region there is a one- 
to-one correspondence between the power-like dependence on this cutoff and the 
position of IR renormalons, which establishes the link between the OPE and the 
structure of the singularities in the Borel plane. However, in the physical region 
we did not succeed to construct a cutoff regularization with this property. In 
other words, for the regularization schemes that we have investigated the cutoff 
dependence does not match with the position of IR renormalons. Although we 
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cannot prove at present that no better regularization exists, we cannot exclude 
the possibility that the ambiguities encountered in our analysis of R e + e - and 
R t hinder the standard application of the OPE in the physical region. This 
would imply that to the resummed perturbative series one would have to add 
nonperturbative power corrections that are not related to vacuum expectation 
values of local, gauge-invariant operators. The values of these nonperturbative 
terms would depend on the resummation prescription chosen for the perturbative 
series. It is not obvious to us which theoretical framework could be employed to 
accommodate these terms. 

It is of great importance to settle the question if this pessimistic possibility 
can be ruled out. Let us illustrate some phenomenological implications with two 
examples. We have mentioned in the introduction that the inclusive B —> X e D e 
decay rate is subject to large two-loop corrections of order /3q a 2 , and that the 
corresponding BLM scale is very low. In the case of b —» u transitions, for 
instance, one finds /U B lm — 0.07 up, in the MS scheme p6|. This observation casts 


doubt on the reliability of the perturbative expansion and calls for an analysis 
of higher-order terms in the series. It is thus very interesting to investigate the 
resummation of renormalon chains in this case. The principle value of the Borel 
integral corresponding to the perturbative series for the inclusive semileptonic 
decay rate has recently been calculated by Ball et al. [|36|. However, based on our 


analysis we expect that also in this case the result of the resummation will not 
be unique, i.e. there will be ambiguities of order 1/m^ or even 1/mb. bi the latter 
case this would limit the accuracy with which one can extract the elements V u b 
and V c b of the Cabibbo-Kobayashi-Maskawa matrix from the analysis of inclusive 
decays of B mesons. 

Our second example concerns the quantity R r . measurements of which have 
been suggested as a way to extract with high precision the strong coupling con- 

m 


stant a s {m 2 ) 


40|j. From the most recent analysis of the exper¬ 


imental data, one has obtained the value R r = 3.638 ± 0.017 and extracted 
a s {ml) = 0.367 ± 0.018 [fnj. The theoretical uncertainty in this result is quoted 
as 5a s (m 2 ) = 0.014 [|33|. To achieve such a high precision, it is essential to rely 


on the standard form of the OPE, in which the magnitude of nonperturbative 
contributions is linked to the structure of renormalon singularities in the Borel 
plane. Using the SVZ approach to relate these nonperturbative contributions 
to vacuum condensates, one finds that the contribution of the gluon condensate, 
which is proportional to (a s G 2 )/m/, is suppressed by two additional powers of 
a s (it is absent in the large-/?o limit considered here), so that the dominant non¬ 
perturbative effects appear at order 1/m® |24|]. Several authors have questioned 
this assumption and investigated, with different conclusions, whether R r could 
receive larger nonperturbative corrections. Both IR and UV renormalons have re¬ 
ceived main attention in these studies [T]J, [42//46]. In the present work we have 


identified a new kind of nonperturbative effects, which are not related to vacuum 
condensates. We have seen that there exist different ways to resum renormalon 
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chain contributions. Some of these resummations are disfavoured for physical 
reasons, but others cannot be distinguished on physical grounds. In the case of 
R T1 we see no physical argument to prefer one of the two resummations T™ s cle and 
Teoreij which differ by terms of order 1/m^. More generally, if the Borel integral 
is regulated by introducing an IR cutoff, one is forced to add a nonperturbative 
contribution of order 1/ml to cancel the cutoff dependence of the perturbative 
calculation. 

To illustrate the phenomenological impact of our results, let us compare the 
various predictions for the perturbative contribution to R r . We write 


1? 7 


= 3(| V ud \ 2 + | V us \ 2 ) Sew {l + few T feiert T feionpert ^ j 


(103) 


where Sew — 1.0194 and <5ew — 0.0010 are known electroweak corrections 
|32[| , (5 P ert is the perturbative correction, and (feonpert denotes the sum of all power- 
suppressed corrections including quark-mass effects, which have been estimated 
to be (fecmpert = —0.019 ± 0.007 |24|]. Accepting this number, one obtains <5^ = 


0.21 ± 0.01 from the experimental measurements. It is conventional to consider 
the perturbative correction as a function of the coupling constant a s (m %) defined 
in the MS scheme. We use the known exact values of the coefficients up to order 
af J22], |23j to correct our result for contributions which are subleading in the 
large-/fe limit and thus not included in our approach. This leads to 


^pert[a S fe«r)] = P[a s ( m r)] “ 1 + 0.0831 


a < 


K) 


2.4115 


7T 


a Amt 


7r 


(104) 


where for T we may substitute one of the resummations T™ cle or Teorei, which 
depend on a s (ml ) through the QCD scale parameter Ay. 

In Fig. H we compare our results with the resummation of Pich and Le Diberder 
as well as with the exact order-a^ expression. There are significant differ- 
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ences between the various predictions. For instance, with the central “experimen¬ 
tal” value 5p ert = 0.21 one would obtain a s {ml ) ~ 0.30 using the principle value 
resummation, a s (m 2 r ) ~ 0.34 using our resummation T™ cle , and a s (m 2 T ) ~ 0.37 
using the resummation of Le Diberder and Pich. If 5p ert is larger, the differences 
even increase. We conclude that the theoretical uncertainty in the perturbative 
contribution to R r is probably larger than estimated in |33|. A reasonable value 
seems to be 

Sa s (ml ) ~ 0.05. (105) 
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Figure 8: Theoretical predictions for the quantity <5 per t: this work, 

eq. ( |68| ) (solid line); principle value of the Borel integral, from || (dash- 
dotted line); resummation of Le Diberder and Pich, from |33|] (dashed 
line); exact order-ct^ result (dotted line). 


Appendix 


Beneke and Braun have developed a different way to resum renormalon chains 
[]H (see also ||). They start from the integral representation (|78| ) for the Borel 
transform and show that the Borel integral performed above the real w-axis is 
equivalent to 


1 1 

S B ovei(M 2 ) = 1 + J dr W(t) Fi(a{rM 2 )) + — JdrW(T-ie), (106) 


1 

A) 


-TL 


where F\ (a) has been dehned in (|34|) . In this approach, the renormalon ambiguity 
is given by 

1 } 

—— 1m / drhF(r — ie) . (107) 

7rp 0 J 


A = 


The function W{r) can be calculated by performing a one-loop calculation 
with a finite gluon mass m 2 g = tM 2 . Its relation to the distribution function 
w{t ) of our approach has been given in (|5T]) . Using this relation, one can recover 
our linear representation (]T|) from the non-linear representation ( |106| ) of Beneke 
and Braun. Inserting ( ^I| ) into ( |106| ) and changing the order of integration (which 
is allowed since the integrals converge), one finds that 


1 

S B ove\{M 2 ) = 1 + — J d tw(t) H(a(rM 2 )), 


(108) 
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where 


H(a) = 




(1 + z) 2 7r 


arccot 


In z + 1/a 


7T 


+ 


3 — 1 / a 


1. 


ie 


(109) 


After an integration by parts in the first term this can be written in the form 

+oo 


H(a) = 


1 - e~ x ! a - 


da; 


ie 


(. x 2 + 7t 2 )(1 + e x l / a ) 


( 110 ) 


The x-integral can be performed using Chauchy’s theorem. There are simple 
poles at x = ±77r and x = 1/a + (2n + 1) nr. We find 
+°° .. . ™ 


da: 


(x 2 + 7T 2 ) (1 + e x ~ 1 ! a ) 1 


g-i/a 

1 


1 — e 1 / a 


E 

n =0 

+ a, 


which leads to 


H(a ) = a + nr 5(l/a) = 


l/a + 2nz7r 1/a + 2(n + 1) in 

(in) 
( 112 ) 


1 /a —ie lnr — Inry — ie 

This shows the equivalence of the two methods in the euclidean region. 

Using the relation ([ST]) , let us derive some properties of the function W(t). 
We note that W(t) is analytic in the complex f-plane with a branch cut along 
the negative axis. For the discontinuity across the cut, we obtain 


- W(r + ie) — W(t — ie) = — ©(—r) — tw(—t) 

2m 1 J dr L 


(113) 


Using this result, we find that our definition of the renormalon ambiguity in (|T1|) 
is in fact equivalent to (|107|), since 


A,S ren = 


tt(3 0 


Im 


U 

J dr W(t - 'V - Z- w[j L ) 


(114) 


For the cases of R e + e - (s) and R T , the functions W e + e - (r) and W T {r) have been 
given in ([16]) and (|57|) . These results agree with the corresponding expressions 
derived in Ref. [|J. For completeness, let us also quote the result for the function 
Wd( t) in the case of the U-fimction. Starting from relation (|8l|) and using the 
result for the distribution function Woij) given in (2l|), we obtain 


W D (r) = 8 C f 


7 5 r , 1 + r 

-In-h 

2 2 r 

r 111 x , 
da;-In 


(^E_2^j ln(l + t) + 2 L 2 {-t) + tL 2 (-1/t) 


- (1 - 2r) 


2 

x(x + r) 

1 + X (1 + XT ) 

17 3 


= 8 C F 14 - 3C(3) + ^6C(3) - — + - lnr) r + 0(r 2 )|. (115) 
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The first two terms in the expansion of this function have been calculated in 
Ref. ||; the exact expression is new. 
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